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We consider a pair of atoms in an arbitrary trapping potential in the presence of magnetically 
tunable Feshbach resonance. We find the energy levels and occupation of the bound molecular states 
taking into account possible coupling between center of mass and relative motion induced by the 
trap. As a specific example we discuss the case of different atomic species in harmonic potential, 
where each atom feels different trapping frequency. 
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I. INTRODUCTION 

Degenerate quantum gases manipulated by electro- 
magnetic fields provide an opportunity to perform quan- 
tum simulations of condensed matter models as well as 
quantum information processing [TH3]. In experiments, 
it is possible not only to precisely control the trapping 
potential, but also to vary the interactions strength, de- 
scribed by the scattering length. This can be done us- 
ing magnetically tunable Feshbach resonances, which are 
one of the most important phenomena in this field ^ [S] . 
Feshbach resonances are a result of coupling between the 
free atom pair and a bound state in the closed chan- 
nel. The position of this state can be controlled using 
an external magnetic field, causing a resonance in the 
scattering length when it crosses the threshold. A large 
scattering length greatly increases three-body losses in 
the system, which was the first sign of such resonances 
in experiment [6 . Nowadays, Feshbach resonances are a 
crucial tool to produce ultracold molecules, which can be 
done e.g. by using time-dependent magnetic fields [THS]. 
Feshbach resonances allow also to observe atom-molecule 
coherence |10] . BEC-BCS crossover pn fT2] and to sim- 
ulate identical fermions with a pair of distinguishable 
atoms [13]. 

The simplest description of Feshbach resonances uses 
a single-channel model in which the atoms are as- 
sumed to interact via a pseudopotential with scatter- 
ing length a depending on the magnetic field [1]. This 
approach is limited to the so-called open-channel domi- 
nated resonances |4j . To increase the precision, especially 
for closed-channel dominated (also called narrow) reso- 
nances, more complex methods are needed. Multichan- 
nel calculations and experimental Feshbach spectroscopy 
have been performed for numerous cases ( |14H18) and 
many more). Full coupled channel calculations can be 
simplified to effective two-channel models [TWET] . 

In the presence of external harmonic traps, single- 
channel pseudopotential calculations can be performed 
analytically [22], but their validity is limited to the case 
when the range of the potential is much smaller than 
the trap width and the scattering length a is sufficiently 
small, so that ka <^ 1 ^3\. Close to the resonance, 
when a can be arbitrarily large, the description has to be 



extended by introducing an energy-dependent scattering 
length [21M1]. 

In free space and in a harmonic trap the center of mass 
and relative motion of a pair of atoms are decoupled, 
and the Feshbach resonance changes only the properties 
of the relative motion. However, for a variety of traps 
used in experiments, such as optical lattices or double 
wells, such separation is not possible anymore. This leads 
to novel phenomena, such as anharmonic confinement- 
induced resonances [27H29] , and formation of states with 
nontrivial angular momentum correlations |30j . Anhar- 
monic terms influence the energy of the pair of atoms and 
the structure of bound states , also in a waveguide [52] 
and in optical lattices [33 37 . 

In this work we consider a pair of atoms confined in 
an arbitrary external trapping potential in the vicinity 
of a Feshbach resonance. We describe the resonance us- 
ing a two-channel configuration interaction (CI) model. 
In our approach we treat the closed channel molecular 
state as a pointlike particle. This approximation results 
in divergencies, which can however be renormalized, in 
close analogy to the free space problem [20]- A conve- 
nient way to do this is to introduce the renormalized 
resonance shift [SH [37] . We then analyze the Feshbach 
resonance in harmonic trap. In this case the renormal- 
ization procedure is particularly simple and can be im- 
plemented numerically without difficulties. This result 
makes it possible to perform efficient numerical calcula- 
tions for the large class of traps where harmonic oscil- 
lator solutions can be used. As a simple example of a 
nonseparable problem, we describe association of a het- 
eronuclear molecule. The separation of center of mass 
and relative motion does not occur, because atoms with 
different masses and polarizabilities feel different trap- 
ping frequencies [301 131j . Because in experiments involv- 
ing mixed species the atoms are always in external traps 
usually having different frequencies, taking the effects of 
nonseparability of the trap into account is crucial for the 
accuracy of calculations. 

Our paper is organized as follows. In Section II we 
outline the basic physics of Feshbach resonances and the 
parameters which describe them. In Section III we gen- 
eralise the theory to the case of arbitrary trapping poten- 
tials and obtain self-consistent equations for the energy 
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levels. Section IV is dedicated to describing Feshbach 
resonances in isotropic harmonic traps and discusses the 
applicability of the method to more complicated cases. 
In Section V we discuss application of our formalism to 
the problem of two different atoms in a harmonic trap. 



II. TWO CHANNEL MODEL OF A FESHBACH 
RESONANCE 

In our description of Feshbach resonances we follow the 
two-channel CI model [71 18] . In this section we briefly re- 
view its characteristics in free space. Let us consider two 
atoms of mass mi and m2- The Hamiltonian consists 
of the open collision channel describing a pair of atoms 
in the spin state |x), the resonant molecular state \n) in 
the closed channel and an interchannel coupling, which 
depends only on the distance between the atoms. By 
introducing the center of mass and relative motion coor- 
dinates 



R 



miTi + m2r2 
mi + 1712 
r = ri - r2, 



(1) 
(2) 



in the absence of external potential the problem can be 
separated. The center of mass solutions are just the plain 
waves. The Hamiltonian of the relative motion reads 



H = |X> (Xl 



2^ 



+ \n) {n\ 



2^1 

(|X> {n\ 



7^ + Utg{r) 

- Urnol{B,r) 



(3) 



Here /i is the reduced mass, U is the background potential 
between the atoms away from the resonance and Wnx is 
the coupling. The CI wave function is given by 

B, r)) = \x) C(e, B)^,{v) + \n) S)$„,o/(r), 

(4) 

where A and C are the amplitudes and <^{r) are 
the channel wave functions, which obey single channel 
Schrodinger equations: 



-^V^ + UM]^,{r)^e^,{r), (5) 



2/i 



U,noi{B,r) ^,noi{r) - v{B)^^oi{r). (6) 



The resonance is controlled by an external magnetic field 
B. v{B) is the energy of the molecule shifted by the 
presence of the field. The effect of the magnetic field 
is not so much to change Umoi, but mostly to control v. 
Close to the resonance, v may be expanded to first order, 
giving 



where s is the difference of magnetic moments between 
the open and closed channel states and Bq is the value 
of the magnetic field at which the energy of the closed 
channel crosses the dissociation threshold in the open 
channel. Other important parameters are the resonance 
width A and background scattering length a^g , connected 
by the identity [3] 



sA 



2kabg 



(8) 



where F is the decay width, given by F(e) = 
27r I {^moi I \ ^<^)\ ■ Within the single-channel descrip- 
tion, the scattering length can be obtained, given by the 
well-known formula [5] 



(9) 



a{B) = flbg 1 



B — Br 



This effective expression for the scattering length does 
not contain the parameter s. 



HI. FESHBACH RESONANCE IN A TRAP 

If the system is in an external trap, the above descrip- 
tion needs to be adjusted. First of all, the separation of 
center of mass and relative motion may no longer be pos- 
sible. We thus rewrite the full Hamiltonian, adding the 
trapping potential Utrap to the interaction U and Umoi- 
The Hamiltonian takes the form 

H = |x> (xl (T + Ubg{r) + Utrap{ri,r2)) + 
+ \n) {n\ (T + U„,,i{r) + {/t™p(ri, ra)) -I- (10) 
+ (|X) {n\ + \n) {x\)Wnx{rl 

where T — "jmT^i ~ 2m^^2 is the kinetic energy oper- 
ator. The general wave function is given by 

r)) = Ix) C^<^^i^^ + |n) Ak^k.^ol{R-, r), 

i k 

(11) 

where the channel wave functions obey 



;j2 f^2 \ 

'^^r~ l^^l + Ubg + Utrap] $^ (r , R) = 6, (r , R) 

^ (12) 



2M 



U, 



rnol 



trap I mol 



;(r,R) 



v{B)^s{B-Bo), 



(7) 



= (i^(S) +efe)$fe„oi(r,R). 

(13) 

Previously it was sufficient to consider only one molecular 
state, but in the presence of the trap one may expect cou- 
pling between different free and molecular levels induced 
by the trapping potential. It is convenient to separate 
the resonant energy shift from the energy of trap exci- 
tations. We will assume that the molecule is a pointlike 
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particle of mass M = mi + m2, which is a reasonable 
assumption as long as the interatomic distance is much 
smaller than characteristic trap lengths. Then the molec- 
ular wave function $fc„jo/(R, r) can be replaced by its 
value at r = and a Dirac delta in r. The i?-dependent 
part satisfies 



2M 

(14) 

Here Utrap(R) = Utrap{r = 0,R). Applying the 
Schrodinger equation H lip) — E lip) to this problem gives 

EC,^e,C, + Y,Vi*Ai, (15) 
I 

EAk = {v{B) + ek)Ak + ^ 14jC,, (16) 



where Vki 



k mol 



The collision takes 



place at short range in comparison to characteristic trap 
lengths, which justifies using only one coupling Wn^- By 
substituting (15) into (161 we get a self-consistent for- 
mula for E: 



{E- 



(17) 



This sum may be divergent, because we treated the 
molecular state as a pointlike particle with 5{r) in the 
relative coordinate, v needs then to be renormalized. 



IV. ISOTROPIC HARMONIC TRAP 

In this section we apply our formalism to the simplest 
possible case of two atoms in an isotropic harmonic trap. 
The trapping potential ^mui'^rf separates center of mass 
and relative motion. We neglect the background interac- 
tion in the open channel, assuming that the background 
scattering length is small The i?-dependent part of 
the problem disappears from the equations, meaning that 
the resonance will not affect the center of mass motion of 
the pair. Furthermore, due to the form we assumed for 
the molecular wave function, only states with £ ~ will 
couple to the closed channel. We thus have 

= IX) E^^^ W + I") ^'^rr^ol{r), (18) 



where 



(dho) \/ Tr^+^j2) normalization factor and 



r(i+3/2) 

^moi{r) is approximated by 5^^^ (r). Additionally the pair 
is described by some center of mass wavefunction which 
does not contribute to the resonance properties. 

The coupling between open and closed channel in a 
harmonic trap may be calculated analytically, using the 



2 r(3+3/2) ^, 



property L^^^(O) 

V, = I d'rMr)Wir)S^'\r) «^/^^|^^. (19) 

where a is a constant. The method to calculate its value 
in terms of experimentally accessible parameters is given 
in the Appendix. By inserting this into ( 17 ) and denoting 
x^{E- 3hu}/2)/2Huj, we get 



E-v 



-a 



E 



r(n-K3/2) 1 
T{n + 1) n — X ' 



(20) 



where the energy of the molecular state was inserted into 
V. The sum in Eq. (20) is divergent. In the numerical 
calculations, when one uses a finite basis, this results in 
dependence of the resonance position on the basis size 
n*. This can be avoided by renormalizing the parame- 
ter V. The divergence can be extracted by adding and 
substracting \l \/n + 1 under the sum. It can be shown 
that 



E 



1 



C(l/2) 2Vn* 



(21) 



where C is the Riemann zeta function. By introducing 
'r(n + 3/2) 1 1 



W{x) = E 



n=0 



r(n +1) n — X + 1 



C(l/2), 



(22) 
(23) 

We may now introduce the renormalized resonance shift 
V* — V ~V 2o?\fri^ which makes the equations convergent 
and ensures that the basis size will not affect the reso- 



which is convergent, we obtain 

E-v^ -o?W(x) - 2a'^Vn*. 



nance properties. In this form Eq. ( 23 ) is very convenient 



for numerical calculations. One can also perform the cal- 
culations without going to the center of mass frame and 
using cartesian coordinates, obtaining the same renor- 
malization condition, as done in [36 in the context of an 
optical lattice. We note that lower dimensional problems 
can also be treated with similar approach .'i!) . 



A. Application to nonseparable problems 

The results obtained in the previous paragraphs are 
not limited to the pure harmonic potential. Instead, they 
can be used to solve a wider class of problems. Let us 
now consider a general trapping potential. The eigen- 
functions in both channels can be expanded in the basis 
of harmonic oscillator states 

$fc(R, r) = E a'l(l)NLM{'R)(l)nim{r) (24) 

i 

<i>r'(R) = Es''/'A^'i'M'(R), (25) 
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where the index i in the sums denotes summation over all 
possible states \N LMnlm) of the pair and j over molec- 
ular states \N'L'M'). We note that only terms with 
L = £ = will couple to the resonance. By using this 
basis to perform the calculations, Eq. ( 17) takes the form 



k' 



N,N',n,n't 



''Nn"'N'n"^N^N 



„k' „k* ^n' Vn 

E - et 



(26) 

where et (e/c) denotes the eigenenergies of the open 
(closed) channel and only expansion coefficients a, c with 
zero angular momentum are present in this formula. As 
the form of the coupling is the same as in the pure har- 
monic oscillaotr case, this equation can be renormalized 
in the same way. This method will be particularily useful 
when one of the following conditions is met: 

• the trapping potential can be described by har- 
monic term plus some perturbation with finite 
strength and range; then the high energy eigen- 
states will not be affected by the perturbation, or 

• the coupling between center of mass and relative 
motion mixes only the states lying close to each 
other. 




FIG. 1: Energy levels in units of the relative trapping fre- 
quency ij for a K-Rb Feshbach resonance. The ratio of the 
trapping frequencies between K and Rb atoms is assumed to 
be 1.4 as in 131]. Five molecular bound states are taken into 
account. The inset shows a closer view of the region where 
the bound states cross the trap levels. 



In both cases a reasonably small basis can be used for 
numerical calculations. 



V. EXAMPLARY APPLICATIONS 

A. Two different atoms in harmonic trap 

As a simple example of a system where the center of 
mass and relative degrees of freedom are coupled, we con- 
sider a combination of two different species with masses 
mi and m2 in a harmonic trap. Due to different masses 
and polarizabilities of the atoms, each atom feels differ- 
ent trapping frequencies wi and uj2. The Hamiltonian of 
the open channel reads 



where fl = 



1 



(27) 



(miujj+mgw^) 
M ' ^ 



and the 



coupling term C in the Hamiltonian is given by [501 EI] 

C = ^l{Lol-uJ|). (28) 

Due to rotational invariance of the Hamiltonian, the total 
angular momentum J of the pair is conserved and we may 
choose it to be equal to zero. The open channel wave 
function may then be expanded in the basis of J = 
harmonic oscillator states I5D1 



i/'Ar^„(R,r) = 



(-1) 



■m——i 



$7Vf™(R)(/)„f(_„)(r), (29) 



where <&(R) is the center of mass harmonic oscillator 
wave function, 0(r) is the relative motion wave function 
and we used the fact th|Lt_the Clebsch-Gordan coefficients 
{liM i.2'm\0Q) give ''"2^1+1 ^f-ii-i- Matrix elements of the 
Hamiltonian ( 29 1 in this basis can be computed analyt- 
ically [30 (only the last term in (29) is not diagonal in 
this basis) . The closed-channel wave function is a super- 
position of ^ = eigenstates <!>„ of a harmonic oscillator 
with frequency Vl and mass M . We thus have 

^/(R, r)) = Ix) CfeVfe(R, r) + \n) ^ Afe$fe(R)5(r), 

k 

(30) 

V'Ar£n(R, r) are the eigenstates of 
(|29). Only the £ = Q components 



k 



where Vfc = Y.Ntn ^Nly^ 

the full Hamiltonian 
are coupled with the closed channel and the coupling has 
the same form as in the previous case of single harmonic 
oscillator, so we can now substitute the wave functions 
and couplings into Eq. ( |17[ ) and solve it numerically. We 
now analyze some particular examples of heteronuclear 
Feshbach resonances. Fig. [T] shows the energy levels in 
the case of a Feshbach resonance between K and Rb at 
547 G [31) . assuming the trapping frequency for rubid- 
ium wi = 10 kHz and for potassium uj2 = 14 kHz. The 
parameters of the resonance can be found in 4J. Away 
from resonance the eigenstates do not contain any bound 
levels. The energies for this case are shown on Fig. [2] 
We note that due to coupling of motional degrees of free- 
dom induced by nonzero C, the eigenstates are composed 
out of several harmonic oscillator levels with different an- 
gular momenta and the corrections to the eigenenergies 
with respect to the uncoupled case are significant. Close 
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FIG. 2: (color online) Energy levels away from the K-Rb res- 
onance. Dashed blue lines depict the uncoupled C = case, 
where the eigenstates can be labeled by the quantum numbers 
Nln. The actual eigenstates (black solid lines) are composed 
from them. 



to the resonance one can see the deeply bound states to 
which one can assign the quantum number N labelling 
the trap level. Then the molecular bound states cross 
with the free atomic states, as shown by the inset of 
Fig. [l] Due to the different symmetry of the states, we 
can expect that these are true level crossings. To verify 
this, we checked numerically that the crossing states are 
orthogonal. 

In Fig. |3] we present the case of Li-Cs resonance at 
B = 816 G, which has recently been observed experi- 
mentally [T5]. Here we assumed that the trapping fre- 
quency is oji = 1 kHz for Cs atoms and ljJ2 — 1.8 kHz 
for lithium. Due to the large mass difference, the ratio 
of trapping frequencies here is bigger than in the K-Rb 
case. As a resul, in the former case the energy levels away 
from the resonance tended to form groups, but here it is 
not the case. Instead we get an energy spectrum which 
looks more complicated, but has similar nature as before. 
The coupling of center of mass and relative motion oc- 
curs to have less impact than in the K-Rb case, so away 
from resonance the eigenstates are less distorted from the 
pure il^Nin states (see Fig.[4|. 



VI. CONCLUSION 

In this paper, we presented a general formalism de- 
scribing Feshbach resonances in an external trap. Our 
method works for both open- and closed channel dom- 
inated resonances and can be applied to nonseparable 
traps. It is particularly efficient when the single par- 
ticle trap eigenf unctions can be expanded in harmonic 
oscillator basis. We presented results for the calculation 
of energy levels of a pair of different atoms, where the 
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FIG. 3: Same as on Fig. [l] but for Li-Cs resonance, where the 
ratio of trapping frequencies is assumed to be 1.8. 
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FIG. 4: (color online) Same as on Fig. [2j but for Li-Cs res- 
onance. The coupling term plays less important role than in 
the K-Rb case. 



trapping frequencies cannot be assumed to be the same. 
Apart from static cases, our formalism allows for calcu- 
lation of the dynamics of the wave functions where the 
trap parameters or magnetic field are changing in time. 
This can be useful for example for quantum computa- 
tions, where control of the qubits will be enhanced by 
Feshbach resonances. 

This work was supported by the Foundation for Pol- 
ish Science International PhD Projects Programme co- 
financed by the EU European Regional Development 
Fund, AQUTE, SFB/TRR21 and National Center for 
Science Grant No. DEC-2011/01/B/ST2/02030. 
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Appendix A: Relation between coupling constant a 
and resonance parameters 

The connection between a and experimentally acces- 
sible parameters can be found by comparison of Eq. |23] 
with the energy of a weakly bound molecule in free space, 
where E = and a{E) is the effective energy depen- 

dent scattering length [4]. This can be done by taking 
the limit a; — >■ at constant E. In this limit x — >■ — cxd 



and W{x) — > — tt-^/— .t. Equation (23) takes the form 



E 



3/2 



-E. 



(Al) 



here was rescaled as <5^a^^, where Uho — Ex- 
panding a{E) into power series according to the effective 
range theory [4] 



1 _ 1 

a{E) abg 



1 2^E 



(A2) 



where tq is the effective range parameter and a^g is the 
s-wave scattering length away from the resonance, we 
obtain the equation 



2^ 



1 



1 



2fi abg 



J2h 



roE. (A3) 



Comparing the energy-dependent and independent 
terms, we conclude that 



2h^ 



1 



(A4) 
(A5) 



We notice that the effective range for this problem is 
negative. It can be found using the definition of energy- 
dependent scattering length in the presence of Feshbach 
resonance [211 125] 



where E^ is the bound state energy. At i? « this re- 
duces to the common formula a — a^g ^1 — b-Bq ) ' 
panding (A6I in E, we get 



1 EbA + {B - Bo)sA - sA^ 



a{E) abg abgEbs{B - Bq - A) 



E 



(A7) 



Comparing this with ( A2 1 and assuming that B is close 
to Bq, we get 



^iabgs{B - Bo) 



(A8) 



Inserting this result into (|A4| and ( A5 1 and neglecting 



the contribution from background scattering length (as 
we are close to resonance), we get 



s{B - Bo) 



a = 2h\ 



'obgSA 
/iTT 



(A9) 
(AlO) 



[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[2] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, 

A. Sen, and U. Sen, Advances in Physics 56, 243 (2007). 
[3] J. J. Garca-RipoU, P. ZoUer, and J. I. Cirac, Journal of 

Physics B: Atomic, Molecular and Optical Physics 38, 

S567 (2005). 

[4] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 

Mod. Phys. 82, 1225 (2010). 
[5] E. Timmermans, P. Tommasini, M. Hussein, and A. Ker- 

man, Phys. Rep. 315, 199 (1999). 
[6] S. Inouye, M. R. Andrews, J. Stenger, H. J. Miesner, 

D. M. Stamper-Kurn, and W. Ketterle, Nature 392, 151 

(1998). 

[7] F. H. Mies, E. Tiesinga, and P. S. Julienne, Phys. Rev. 

A 61, 022721 (2000). 
[8] P. S. Julienne, E. Tiesinga, and T. Kohler, Journal of 

Modern Optics 51, 1787 (2004). 
[9] T. Kohler, K. Goral, and P. S. Julienne, Rev. Mod. Phys. 

78, 1311 (2006). 
[10] E. Donley, N. Claussen, S. Thompson, and C. Wieman, 

Nature 417, 529 (2002). 



[11] T. Bourdel, L. Khaykovich, J. CubizoUes, J. Zhang, 
F. Chevy, M. Teichmann, L. Tarruell, S. J. J. M. F. 
Kokkelmans, and C. Salomon, Phys. Rev. Lett. 93, 
050401 (2004). 

[12] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, 

C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. 

Lett. 92, 203201 (2004). 
[13] G. Ziirn, F. Serwane, T. Lompe, A. N. Wenz, M. G. Ries, 

J. E. Bohn, and S. Jochim, Phys. Rev. Lett. 108, 075303 

(2012). 

[14] P. J. Leo, C. J. Williams, and P. S. Julienne, Phys. Rev. 

Lett. 85, 2721 (2000). 
[15] C. Marzok, B. Deh, C. Zimmermann, P. W. Courteille, 

E. Tiemann, Y. V. Vanne, and A. Saenz, Phys. Rev. A 

79, 012717 (2009). 
[16] T. Schuster, R. Scelle, A. Trautmann, S. Knoop, M. K. 

Oberthaler, M. M. Haverhals, M. R. Goosen, S. J. J. 

M. F. Kokkelmans, and E. Tiemann, Phys. Rev. A 85, 

042721 (2012). 

[17] T. Takekoshi, M. Debatin, R. Rameshan, F. Ferlaino, 
R. Grimm, H.-C. Nagerl, C. R. Le Sueur, J. M. Hutson, 
P. S. Julienne, S. Kotochigova, et al., Phys. Rev. A 85, 



7 



032506 (2012). 

[18] M. Repp, R. Pires, J. Ulmanis, R. Heck, E. D. Kuhnle, 
M. WeidemuUer, and E. Tiemann, Phys. Rev. A 87, 
010701 (2013). 

[19] F. H. Mies and M. Raoult, Phys. Rev. A 62, 012708 
(2000). 

[20] S. J. J. M. F. Kokkelmans, J. N. Milstein, M. L. Chiofalo, 
R. Walsor, and M. J. Holland, Phys. Rev. A 65, 053617 

(2002). 

[21] P.-I. Schneider and A. Saenz, Phys. Rev. A 80, 061401 
(2009). 

[22] T. Bush, B. Englert, K. Rzazewski, and M. Wilkens, 
Found. Phys. 28, 549 (1998). 

[23] E. Ticsinga, C. J. WiUiams, F. H. Mies, and P. S. Juli- 
enne, Phys. Rev. A 61, 063416 (2000). 

[24] E. L. Boida, E. Tiesinga, and P. S. Julienne, Phys. Rev. 
A 66, 013403 (2002). 

[25] E. L. Bolda, E. Tiesinga, and P. S. Julienne, Phys. Rev. 
A 68, 032702 (2003). 

[26] Z. Idziaszek and T. Calarco, Phys. Rev. A 74, 022712 
(2006). 

[27] E. L. Bolda, E. Tiesinga, and P. S. Julienne, Phys. Rev. 

A 71, 033404 (2005). 
[28] J. P. Kestner and L.-M. Duan, New Journal of Physics 

12, 053016 (2010). 



[29] S.-G. Peng, H. Hu, X.-J. Liu, and P. D. Drummond, 

Phys. Rev. A 84, 043619 (2011). 
[30] J. F. Bertelsen and K. M0lmer, Phys. Rev. A 76, 043615 

(2007). 

[31] F. Deuretzbacher, K. Plassmeier, D. Pfannkuche, 
F. Werner, C. Ospelkaus, S. Ospelkaus, K. Sengstock, 
and K. Bongs, Phys. Rev. A 77, 032726 (2008). 

[32] V. S. Melezhik and P. Schmelcher, New Journal of 
Physics 11, 073031 (2009). 

[33] D. B. M. Dickerschcid, U. Al Kliawaja, D. van Oosten, 
and H. T. C. Stoof, Phys. Rev. A 71, 043604 (2005). 

[34] H. P. Biichler, Phys. Rev. Lett. 104, 090402 (2010). 

[35] S. Sala, P.-I. Schneider, and A. Saenz, Phys. Rev. Lett. 
109, 073201 (2012). 

[36] R. B. Diener and T.-L. Ho, Phys. Rev. Lett. 96, 010402 

(2006) . 

[37] M. L. Wall and L. D. Carr, Phys. Rev. A 87, 033601 

(2013). 

[38] B. Marcelis, E. G. M. van Kempen, B. J. Verhaar, and 
S. J. J. M. F. Kokkelmans, Phys. Rev. A 70, 012701 

(2004). 

[39] J. P. Kestner and L.-M. Duan, Phys. Rev. A 76, 063610 

(2007) . 



